The Insulating State of Matter: A Geometrical Theory 
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In 1964 W. Kohn published the milestone paper "Theory of the insulating state", according to 
which insulators and metals differ in their ground state. Even before the system is excited by any 
probe, a different organization of the electrons is present in the ground state and this is the key 
feature discriminating between insulators and metals. However, the theory of the insulating state 
remained somewhat incomplete until the late 1990s; this review addresses the recent developments. 
The many-body ground wavefunction of any insulator is characterized by means of geometrical 
concepts (Berry phase, connection, curvature, Chern number, quantum metric). Among them, it is 
the quantum metric which sharply characterizes the insulating state of matter. The theory deals 
on a common ground with several kinds of insulators: band insulators, Mott insulators, Anderson 
insulators, quantum Hall insulators, Chern and topological insulators 
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INTRODUCTION 

The standard textbook approach to the insulating 
state of matter is based on band theory. Following 
Bloch's theorem [1] in 1928, the main result is due to 
Wilson in 1931 [2]. The single-particle spectrum of a 
lattice-periodical Hamiltonian is in general gapped, and 
the electron count determines where the Fermi level lies. 
If it crosses a band one has a conductor: an applied 
electric field induces free acceleration of the electrons (at 
T = in absence of dissipation). If the Fermi level lies 
instead in a gap, one has an insulator: in presence of 
a field the electronic system polarizes, but no steady- 
state current flows for T — > 0. This very successful 
theory explains the insulating/conducting behavior of 
most common materials across the periodic table, for 
which band structure calculations became soon available. 

At the root of band theory are two basic assumptions: 
the electrons are noninteracting (in a mean-field sense), 
and the solid is crystalline. By the early 1960s, however, 
it became clear that there are solids where these two 
assumptions are very far from the truth, and where 
the insulating behaviour is due to completely different 
mechanisms. The works of Mott in 1949 [3] and of 
Anderson in 1958 [4] opened new avenues in condensed 
matter physics. In the materials which we now call Mott 
insulators the insulating behaviour is due to electron 
correlation [5], while in those called Anderson insulators 
it is due to lattice disorder [6]. 

In a milestone paper appeared in 1964 Kohn [7] 
provided a more comprehensive characterization of the 
insulating state of matter, which encompasses band 
insulators, Mott insulators, Anderson insulators, and 
eventually any kind of insulating material. According 
to Kohn, the electrons in the insulating state satisfy 
a many-electron localization condition [8]. This kind 
of localization must be defined in a subtle way given 
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that, for instance, the HamUtonian eigenstates in a band 
insulator are obviously not localized. According to the 
original Kohn's formulation, the insulating behaviour 
arises whenever the ground-state wavefunction of an 
extended system breaks up into a sum of contributions 
which are localized in essentially disconnected regions of 
the many-electron configuration space. 

Kohn's theory remained little visited for many years 
[9] until the 1990s, when a breakthrough occurred 
in electronic structure theory: the modern theory of 
polarization (for historical presentations, see e.g Refs. 
[10-12]). Inspired by the fact that electrical polarization 
discriminates qualitatively between insulators and metals, 
Resta and Sorella [13] in 1999 provided a definition of 
many-electron localization rather different from Kohn's, 
and deeply rooted in the theory of polarization. Their 
program was completed soon after by Souza, Wilkens and 
Martin [14] (hereafter quoted as SWM), thus providing 
the foundations of the modern theory of the insulating 
state. An early review paper appeared in 2002 [15]. 

Several advances occurred afterwards, and the present 
review illustrates most of them. But the major aim of this 
work is to provide a thorough formulation of the modern 
theory of the insulating state by means of geometrical 
concepts, which remained somewhat hidden and implicit 
in the original literature. While the theory of polarization 
is based on a Berry phase, the theory of the insulating 
state as presented here is rooted in quantum distance and 
metric; most of the results previous published in different 
form will be revisited emphasizing their geometrical 
content. 

GEOMETRY IN QUANTUM MECHANICS 

Let us assume that a generic time-independent 
quantum Hamiltonian has a parametric dependence. The 
Schrodinger equation is 

H{K)\^{K))=E{K)\^{n)), (1) 

where the d-dimensional real parameter n is defined in 
a suitable domain of R"^: a 2d k has been chosen for 
display in Fig. 1. In this Section we start discussing 
the most general case, and therefore for the time being 
we do not specify which quantum system is described 
by this Hamiltonian, nor what the physical meaning of 
the parameter k is. Later on, |\1/(k)) will be identified 
with a many-electron wavefunction. The state vectors 
\'^{k)) are all supposed to be normalized and to reside 
in the same Hilbert space: this amounts to saying that 
the wavefunctions are supposed to obey K-independent 
boundary conditions. We focus on the ground state 
|\1/o(/«)), and we assume it to be nondegenerate for k 
in some domain of M''. 

Any quantum mechanical state vector is arbitrary by 
a constant phase factor. Here we refer to choosing this 




FIG. 1: State vectors in the two-dimensional K-space. 
The phase difference between two of them is defined as 
e-'^-^i^ = and their distance as iJfa = 

1-|{*o(ki)|*o(k2))|^ 

phase as to the choice of the gauge. All measurable 
quantities (e.g. expectation values) are obviously gauge- 
invariant, but the reverse is also true: all gauge-invariant 
properties are — at least in principle — measurable. 

It is expedient to define the ground-state projector 
(a.k.a. density matrix) and its complement, i.e. 

P(k) = |vI>o(k))(vI/„(k)|; Q{K) = i-P{K). (2) 

Both P{k) and Q{k) are obviously gauge- invariant. 

Phases and distances 

We define the phase difference between the ground 
eigenstates at two different k. points in the most natural 
way: 

^-iA^,. ^ (^0(<^l)|^0(^C2)) . , . 

|(*0(«l)|*0(K2))| ' ^ ' 

A<^i2 = - Im log (*o(ki)|*o(k2)) . (4) 

For any given choice of the two states, Eqs. (3) and (4) 
provide a A9912 which is unique modulo 2tt, except in the 
very special case that the states arc orthogonal. However, 
it is also clear that such Aipi2 is gauge-dependent and 
cannot have, by itself, any physical meaning. 

The distance between quantum states has been defined 
by Bures [16] as: 

£>?2 = l-K*o(Ki)|*o('*2))r. (5) 

Such distance fulfills the familiar axioms from calculus 
textbooks; it vanishes when the two states physically 
coincide (i.e. independently of the phase factor), and is 
maximum (equal to one) when the states are orthogonal. 
At variance with A(pi2 defined above, the Bures distance 
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FIG. 2: A closed path joining four states in K-space. 



is gauge- invariant, and can be explicitly expressed in 
terms of ground-state projectors 

Dl^ = l- Tr {P{k,)P{k2)}, (6) 

where "Tr" is the trace over the Hilbert space. 

Berry phase 

We have already observed that the phase difference 
A(^i2 between any two states is gauge-dependent and 
cannot have any physical meaning. Matters are quite 
different when we consider the total phase difference 
along a closed path which joins several points in a given 
order, as shown in Fig. 2: 

7 = Aipi2 + Ay>23 + A<^34 -I- A(p4i 

= - Imlog (*o(ki)|*o(k2))(*o(k2)|*o(k3)) X 

X (*0(K3)|*0(K4)>(*0(K4)|*0(Ki)). (7) 

It is now clear that all the gauge-arbitrary phases cancel 
in pairs, such as to make the overall phase 7 a gauge- 
invariant quantity. The above simple-minded algebra 
leads to a result of overwhelming physical importance: in 
fact, a gauge-invariant quantity is potentially a physical 
observable. In essence, this is the revolutionary message 
of Berry's celebrated paper, appeared in 1984 [17, 18]. 

Next we consider a smooth closed curve C in the 
parameter domain, such as in Fig. 3, and we discretize 
it with a set of points on it. Using Eq. (3), we write the 
phase difference between any two contiguous points as 



i(*o(K)i*o(K+AAv))r 



(8) 



If wc further assume that the gauge is so chosen that the 
phase varies in a differentiable way along the path, then 
from Eq. (8) we get to leading order in Ak: 



FIG. 3: A smooth closed curve C in K-spacc, and its 
discretization. 



In the limiting case of a set of points which becomes 
dense on the continuous path, the total phase difference 
7 converges to a circuit integral: 
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A.{n) ■ dK, 



where A{k) is called the Berry connection: 
>l(«) = i(*o(«)|VK*o(«)). 



(10) 



(11) 



Since the state vectors are assumed to be normalized 
at any k, the connection is real; we can therefore 
equivalently write 



^(K) = -Im (*o(k)|V,^*o(/«)>- 



(12) 



Dealing with the disparate manifestations of the Berry 
phase which occur in molecular and condensed matter 
phenomena is clearly beyond the scope of the present 
review; we quote Refs. [18-22] out of many possible 
references. 



Connection and curvature 

The loop integral of the Berry connection (i.e. the 
Berry phase 7) is non trivial in two cases: either the curl 
of A.{k) is nonzero, or the curl is zero but the curve C 
is not in a simply connected domain. In the former case, 
we can invoke Stokes theorem; the formulation is very 
simple when /c is a 3d parameter. If C is the boundary 
of a surface S (i.e. C = 9E), and the curl of A{k) is 
regular on S, then Stokes' theorem reads 



7 = 



® A.{k) ■ dn — r2(/«) • n da, 



(13) 



-iA(p ~ (*o('*)|Vk*o(k)) • Ak. 



(9) 



'as Jn 
where Jl is the Berry curvature, defined as 

"(k) = Vk X A{k) = -Im (V/«^'o(k)1 X ]Vk*o('«)) 
= *(Vk*o(k)| X |Vk*o(k)), (14) 
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with the usual meaning of the cross product between 
three-component bra and ket states. Equation (13) may 
be spelled out by saying that the curvature is the Berry 
phase per unit area of S. 

For d ^ 2> the Berry curvature is conveniently written 
as the dx d antisymmetric matrix 

^^p{K) - -2Im (a„*o(/«)|a^j*o(/«)); (15) 

Greek subscripts are Cartesian coordinates throughout, 
and doi = d/dKa- The Stokes theorem can still be 
applied, generalizing Eq. (13) to 

^ =]- [ dn"' AdK^ n^pin). (16) 

The Berry connection is also known as "gauge 
potential", and the Berry curvature as "gauge field" 
[20]. It is worth pointing out that the former is 
gauge-dependent, while the latter is gauge- invariant 
and therefore corresponds in general to a measurable 
quantity, even before any integration. The two quantities 
play (in K-space) a similar role as the vector potential and 
the magnetic field in elementary magnetostatics: A(r) is 
gauge-dependent, nonmeasurable; B(r) = Vr x A(r) is 
gauge- invariant, measurable. 

The Berry phase 7, defined as the integral over a closed 
curve C of the connection, is gauge invariant only modulo 
2tt. This indeterminacy is resolved by Eqs. (13) and (16) 
whenever the curve C is the boundary of a surface S 
where the curvature is regular. In fact, the curvature is 
gauge-invariant and has no modulo 2-k indeterminacy. 

Chern number 

The rhs of Eqs. (13) and (16) is the flux of the Berry 
curvature on the surface S; such flux remains meaningful 
even on a closed surface (e.g. a sphere or a torus), in 
which case dY, is the empty set. The key result is that 
such an integral is quantized. Here we limit ourselves to 
3d where 

j il{K) ■ -a. da = 2itCx; (17) 

C\ is an integer, called Chern number of the first class. 

The proof is based on a similar algebra as for 
Dirac's theory of the magnetic monopole [19, 23], and 
is also closely related to the Gauss-Bonnet theorem in 
differential geometry. 

The curvature is regular (and divergence-free) on the 
closed surface S; the Ihs of Eq. (17) is the fiux of 
J1(k) across S. The integrand $!(/«) is the curl of the 
connection A.{k)] the latter in general cannot be defined 
as a single- valued function globally on E, but only on 
patches of it [19, 23]. To fix the ideas, suppose that 
r2(K.) is singular at k, = 0, and that S is a spherical 



surface centered at the origin. We cut this surface at the 
equator = and we consider the flux across the two 
open surfaces: 

/ n{K,)-nda= I n{K)-nda+ n{K)-nda. (18) 

We notice that = 9S_ = C, but the surface normals 
n have opposite orientations. From Stokes theorem, 
Eq. (13), we get: 

/ n{K) ■ n da = ± (f> A±{k) ■ dn (19) 
/ fl{K)-nda=f A+{K)-dK-f A-{k) ■ dn . (20) 

JT. Jc Jc 

The two upper and lower Berry connections >t± (/«) may 
only differ by a gauge transformation; the rhs of Eq. (20) 
is the difference of two Berry phases on the same path 
and is necessarily a multiple of 2-k. This concludes the 
proof of Eq. (17). 

We emphasize that the Chern number is a robust 
topological invariant of the wavefunction, and is at the 
origin of observable effects. The Chern number made 
its first appearance in electronic structure in 1982, in 
the famous TKNN paper about the quantum Hall effect 
[24]. Nowadays more complex topological invariants are 
in fashion, and they characterize a completely novel class 
of insulators, called "topological insulators" [25-31]. 

Metric 

Starting from Eq. (5), the infinitesimal distance is 

d 

,K+dK = X/ 9ap{K)dKadK^, (21) 

where the metric tensor is easily shown to be 

QalsiK) = Re (9c«*o(«)|9/3*o(k)) 

- {da^o{K)\^o{l^)){-^oiK)\d0^o{K)) 

= Re {d^^o{ii)\Q{K)\d0^o{K)); (22) 

the projector Q{k) is the same as defined in Eq. (2). This 
quantum metric tensor was first proposed by Provost and 
Vallce in 1980 [32]. 

At this point we may compare Eq. (22) to Eq. (15), 
noticing that the insertion of Q{k) is irrelevant in the 
latter, i.e. 

fl^^in) = -2Im (a„*o(K)10(K)ia^*o(K))- (23) 

It is therefore clear that gajs and ^lajs are, apart for a 
trivial —2 factor, the real (symmetric) and the imaginary 
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(antisymmetric) parts of the same tensor, which we are 
going to call rjafs in the following; 

Va^il^) = R*o(K)|Q(,c)|a^^rQ(K)). (24) 

The curvature, the metric, and hence the full tensor ?7„^ 
are all gauge-invariant quantities. A compact equivalent 
expression is 

77„^(/t) = Tr {dc,PiK)Q{K)dpP{K)}, (25) 

manifestly gauge-invariant and Hermitian. 

Sum over states 

The K-derivatives entering many of the previous 
equations, e.g. Eq. (24), can be expressed starting from 
perturbation theory: 

|vE'o(k + Ak))-|v1/o(k)) (26) 
(*„(«)! [H{k + Ak) - H{k) ] |*o(k)) . 



|a„*o(K)) - E Eo{n)-E4K) • ^'^^ 

Thiese seemingly obvious and innocent formulae need 
some caveat. It is clear that inserting Eq. (27) into 
the Berry connection, Eqs. (11) and (12), we get a 
vanishing resiilt for any k. This happens because 
the simple expression of Eq. (26) corresponds to a 
very specific gauge choice (called the parallel-transport 
gauge [21]); multiplying the rhs by dependent 
phase factor is legitimate, and must not modify any 
physical result, while the Berry connection is instead 
affected. Nonetheless, since our rjapifi) is a gauge- 
invariant quantity, we may safely evaluate it in any 
gauge, including the parallel-transport gauge, implicit in 
Eq. (27). The result is 

ric^iK) (28) 

_ {^o{K)\daH{K)\^4l^)){^nif^)\d0H{K)\^o{K)) 



This expression shows explicitly that both the curvature 
and the metric are ill defined and singular wherever the 
ground state is degenerate with the first excited state. 
Indeed, this is the main reason why the domain may 
happen not to be simply connected. 

MANY-ELECTRON WAVEFUNCTION 

In the previous Section we have discussed some 
geometrical aspects of quantum mechanics, remaining 



at a rather abstract level, and without addressing 
either specific physical systems, or measurable quantities. 
Prom now on, instead, we address a many-electron 
system, whose most general Hamiltonian we write, in the 
Schrodinger representation and in Gaussian units, as 
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(29) 



i=l 



where We and — e are the electron mass and charge, 
respectively (e > 0). Equation (29) is exact in the 
nonrclativistic, infinite-nuclear-mass limit. In Eq. (29) 
A(r) is a vector potential of magnetic origin and the 
potential Y includes one-body and two-body (electron- 
electron) contributions. The 3d parameter k, having the 
dimensions of an inverse length, can be regarded as a 
pure gauge, constant in space, which does not affect the 
fields (electric and magnetic) and therefore is irrelevant 
in classical mechanics. Matters are different in quantum 
mechanics, where the potentials (vector and scalar) are 
the most fundamental quantities, as we know well since 
the Bohm-Aharonov paper in 1959 [33, 34]. In our case 
the parameter k, called "flux" or "twist", affects indeed 
the eigenvalues and eigenvectors depending on the chosen 
boundary conditions. 



Boundary conditions 

Two very different kinds of boundary conditions 

occur in electronic structure theory; "open" boundary 
conditions (OBCs) and periodic boundary conditions 
(PBCs). The former are appropriate to molecular 
physics, and require that the many-electron wavefunction 
of a bound state is square-integrable over the whole 
coordinate space K^^. PBCs are instead appropriate 
for extended systems, either crystalline or disordered, 
either independent-electron or correlated. For the sake 
of simplicity here we adopt PBCs over a cubic box of 
side L, meaning that the eigenstates of Eq. (29) at 
any given k are Born-von-Karman periodic in the cubic 
box over each electron coordinate independently; their 
Cartesian components a arc then equivalent to the 
angles 2TTri^a/L). The potential V enjoys the same 
periodicity, which implies that the microscopic electric 
field averages to zero over the sample. The eigenfunctions 
are normalized in the hypercube of volume L^^ . 

As previously observed, the choice of boundary 
conditions (either OBCs or PBCs) corresponds to 
choosing the Hilbert space, which in turn affects 
profoundly the geometrical properties discussed in Sec. . 



6 



Open boundary conditions 

The case of OBCs is by far the simplest. We write 

for the sake of siniphcity the ground state of Eq. (29) at 
« = as 1*0) = 1*0(0)), and we define the many-body 
position operator as 

N 

r = J2^i. (30) 
1=1 

It is straightforward to verify that the state e~*'* '"|^'o) 
fulfiUs both the Schrodinger equation at the given k 
and the boundary conditions with a K-independent 
eigenvahic. In jargon wc say that the K-depcndencc is 
easily "gauged away" within OBCs; we anticipate that 
matters are quite different within PBCs. 

The state e~''* ''|5'o) coincides therefore with the 
ground eigenstate |\1/o(k.)) of the twisted Hamiltonian, 
Eq. (29). It is legitimate to multiply this eigenvector 
by any K-dcpendcnt (and position-independent) phase 
factor; our choice is then 

|*o(«))=e-"*-(^-'^)|*o), (31) 

where d = (^'o|r|\E'o) is the electronic dipole of the 
molecular system. It follows that the K-derivative needed 
in our geometrical quantities (curvature and metric) is 

|Vk*o) = -i{r - d)|*o) = -iQ{0) r|*o), (32) 

Exploiting the idempotency of the projector, the tensor 
rja/s, evaluated at k = 0, is 

Vap{0) = («'o|^aQ(O)f0|*o) 

= (*o|rar^|*o) - (*o|r«|*o)(*o|r/3|*o)(33) 

This is clearly a real symmetric tensor, hence the Berry 
curvature vanishes within OBCs, and the metric (?a/3(0) 
coincides with r/ap (0) . It will be shown that within PBCs, 
instead, the tensor 77a/3(0) may acquire a nonvanishing 
imaginary part. 

Clearly, the metric tensor in Eq. (33) is the cumulant 
second moment of the position operator, or equivalently 
the ground state fluctuation of the dipole of the molecular 
system; this quantity is extensive (scales like N in 
macroscopically homogenous systems). Wc anticipate 
that r]a/3{0)/N discriminates, in the large A^-limit, 
between insulators and metals. 

Periodic boundary conditions 

First of all, some preliminary observations about the 
position operator are in order. The simple multiplicative 

operator f, as defined in Eq. (30), is "forbidden" within 
PBCs: in fact it maps any periodic wavefunction j^*) into 



the nonperiodic function fl^*). In other words, it maps a 
function which belongs to the Hilbert space to a function 
outside of it [35] . Incidentally, this is the main reason why 
the polarization problem has remained unsolved until the 
early 1990s (for historical presentations, see e.g Refs. [10- 
12]). In the present context, formulae like Eq. (33) are 
ill defined and absurd within PBCs. 

Suppose that j^l/o) is the ground eigenstate of Eq. (29) 
at K = 0. Then the simple state e-''* '|*o) fulfills the 
Schrodinger equation but not — for a general k — PBCs. 
It is therefore not an eigenstate, except for the discrete 
K values 

277 

«mim2m3 = -^(^161 + m2e2 + maBa), (34) 

where TOq, G Z, and Cq, are the Cartesian versors. For 
fractional values of k — i.e. for values different from those 
in Eq. (34) — the /c-dependence of the eigenvalues and 
eigenvectors of Eq. (29) is nontrivial. 

If |4'o(k)) is the genuine ground eigenstate of Eq. (29) 
within PBCs, then the auxiliary function |^o('«)) = 
e*'* '' I *□(/«)) is a solution of H{0), but fulfills quasi- 
periodic boundary conditions: at any two opposite 
faces of the cube the wavefunction differs by a k- 
dependent phase factor. In other words the problem 
can be formulated in two equivalent ways: either 
the Hamiltonian is K-dependent, as in Eq. (29), 
and the boundary conditions arc /^-independent; or 
the Hamiltonian is K.-independent but the boundary 
conditions are "twisted" in a K-dependent way. 

We briefly switch to a Id formulation in order to 
explain the role of k, and even the alternative semantics 
("flux" instead of "twist"). Setting the magnetic vector 
potential to zero, the Id version of Eq. (29) is 

^('^^ = ^ J2^P^ + ^«)' + ^' (35) 
^ j=i 

and we require the wavefunction to be periodical with 

period L over each electron coordinate independently. 
Alternatively, we may regard this problem as if the 
electrons were confined to a circular rail of circumference 
L (sec Fig. 4). There is no fields (electric or magnetic) 
on the rail, but the Hamiltonian is the same as if a 
constant vector potential of intensity A = chu/e were 
present along the rail. This corresponds to a magnetic 
flux $ = chnL/e piercing the surface encircled by the 
rail, in a region not visited by the electronic system; it has 
been appropriately called by some authors "inaccessible 
flux". This flux affects nontrivially the eigenvalues and 
eigenvectors, as it is easily verifled in the trivial case 
V = 0. The inaccessible flux has therefore observable 
effects; a similar algebra in a different context leads 
indeed to the Bohm-Aharonov effect [33, 34]. When k 
assumes one of the values in Eq. (34), then the flux $ 
is an integer multiple of the elementary flux quantum 
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FIG. 4: The electron motion is confined to a circular rail. 

A constant vector potential along the rail, as in Eq. (35), 
corresponds to vanishing fields (electric and magnetic), yet 
the spectrum depends on the "inaccessible flux" piercing the 
surface encircled by the rail. 

$0 = 2iThc/e = hc/e ($o = h/e in SI units). We stress 
that only the fractional part of the flux affects the results 
in a nontrivial way. 

It has already been observed that within PBCs the 
tensor rjap (0) cannot be simplifled — as e.g. in Eq. (33) — 
and must be therefore addressed in its original form, 
Eq. (24), by actually evaluating the k. derivatives. 
As within OBCs, even within PBCs these tensors are 
extensive. In the thermodynamic limit (i.c TV — > oo, 
keeping N/L'^ constant), the well defined quantity is 

Finally, wc observe that for time-reversal symmetric 
systems ?7q,^(0) is real symmetric, and coincides therefore 
with the metric tensor 5q,^(0). 

THE LOCALIZATION TENSOR 

In the milestone paper, appeared in 1964 under the 
title "Theory of the insulating state" [7, 8], W. Kohn 
addressed the insulating character of a solid from a 
very general viewpoint. Kohn focusses on the electronic 
ground state as a whole. His outstanding message is 
that "insulating behavior does not appear to depend on 
the notion of a filled band but rather refiects a certain 
type of organization of the electrons" . Furthermore, the 
electrons "so organize themselves as to satisfy a many 
electron localization condition" [8]. Such localization, 
however, must be defined in a subtle way. For instance, 
the (single-particle) Hamiltonian eigenstates of a band 
insulator are definitely not localized. 

The modern formulations of the theory of the 
insulating state are all based on a localization tensor 
(squared localization length in Id), first introduced by 
Resta and Sorella in 1999 [13], as the main tool to 
address the "organization of the electrons" in the ground 
state. This work was followed soon afterwards by SWM 
[14], where additional results are found, most notably 
relating localization to conductivity. Since then, most 
authors (including the present one) have adopted the 



notation {rarp)c for the localization tensor, where "c" 
stays for cumulant. The formulation within OBCs, using 
the language and the notations familiar in quantum 
chemistry, dates since 2006 [36]. 

The tensor (rar^)c has the dimensions of a squared 
length; it is an intensive quantity that characterizes the 
ground-state many-body wavefunction whole. Its 

key virtue is that it discriminates between insulators and 
metals: it is finite in the former case and divergent (in 
the large-system limit) in the latter. This is the main 
message of the present work (and of the modern theory 
of the insulating state): we are going to prove it below. 
Sect. 

Several definitions of the localization tensor have been 
given in the literature, all of them equivalent. Here 
we give a novel one, deeply rooted in the geometrical 
concepts discussed so far: the localization tensor is the 
intensive quantity 

(rar/3)c -??a^(0)/iV, (36) 

where the thermodynamic limit is understood. A glance 
at the OBCs expression, Eq. (33), explains the reason for 
the notation, which wc adopt within PBCs as well. 

Until 2005 the theory of the insulating state implicitly 
addressed time-reversal invariant systems only, where the 
T] tensor is real symmetric, and coincides with the metric: 
in such systems therefore 

(rar/3)c = 9a0{O)/N. (37) 

It was found in 2005 [37] that in absence of time- 
reversal symmetry and within PBCs — the tensor {rarf))c 
is naturally endowed with an antisymmetric imaginary 
part, whose physical meaning is also outstanding. This 
is discussed in Sects, and . 

Discretized formulae 

While practical implementations within OBCs invari- 
ably use Eq. (33), those within PBCs require a discretiza- 
tion in M-space in order to evaluate the derivatives ap- 
pearing in Eq. (24) . A key feature is that any discretiza- 
tion must be numerically gauge invariant. Here we are 
going to present a novel discretized formula, based on 
Eq. (25) at /« = 0. 

By definition we have 

a„P(0) = hm i[P(fceJ - P(0)], (38) 
hence for small k 

VMO) ^ ^Tr{[P(fce„)-P(0)]g(0)[P(fce^)-P(0)]} 

= -^Tr {P(te„)Q(0)P(fce^)} (39) 
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^^Q-jwhere the rhs is nothing else than h times the velocity 
operator v; Eq. (43) becomes then 



- ^ (*o (A;e„) I *o) (^-0 1 *o (fce^ )) (*o (fce^) | *o (Ae„ )) . 

The gauge invariance is pespicuous: all the arbitrary 
phase factors cancel in pairs. 

Equation (40) is still a three-point formula; it is 
transformed to a single-point one using a key result 
from Ref. [13]. Whenever k is equal to ^ times an 
integer, the wavefunction e~*'^''° |\E'o(0)) obeys PBCs — 
see also Eq. (34) — and is an eigenstate of If(fce„) with 
eigenvalue Eq. For insulating systems the ground state is 
nondegenerate even in the thermodynamic limit, hence 



l*o(^e«)) 



»/^|vl/o). 



(41) 



apart possibly by a phase factor, irrelevant here. 
Replacing k with ^ for large L in Eq. (40), and using 
Eq. (41) we get the sought for formula, where a single 
diagonalization of Eq. (29) at /« = is enough to obtain 
both the real and imaginary part of the localization 
tensor {raTp)c- Notice that the first line of Eq. (40) is 
real. 

The formulae previous reported in the literature are 
inspired by the Berry-phase formulae — see Eq. (7) — and 
based on logarithms. To retrieve them, it is enough to 
use log(l -t- e) ~ e in Eq. (40): 



r, rm-. ^ In:^ (^o(fce„)|^o(fce^)) 



(42) 



which provides (rQ,r^)c identical in form to Eq. (86) in 
Ref. [15]; for a = /3 the formula is also identical to the 
archetypical one, proposed by Resta and Sorella in 1999 
[13]. For the real part of (rQ,r^)c the logarithm formula, 
Eq. (42), has some advantages over Eq. (40). For the 
imaginary part, instead, the second line of Eq. (40) 
provides a much better result, since it is not affected 
by any modulo 2-k indeterminacy. Getting rid of such 
indeterminacy, in fact, is essential when dealing with 
Chern numbers (Sects, and ). 



Sum over states again 



We express (r^r^l 
Eq. (28) : 



using the sum-over-states formula, 



(*o|5aff(0)|*„)(«'„|9/3^?(0)|fo) 



The K-derivative of the Hamiltonian of Eq. (29) is 



V'.^(0) = -E[p^ + ^A(rO 



(43) 



(44) 



{rar0)c = 



(^o|i}a|^„)(^„|f7g|^o) 

{Eo - Er,r 



= -y' 



(*o|^a|«'n)(«'«|v/3|*0> 



(45) 



On 



where won = {En - Eo)/h. 

The velocity operator is also commonly expressed as 
V = i[H(0),r]/h, but it is worth emphasizing that while 
the position r is well defined within OBCs and ill defined 
within PBCs, the velocity v is well defined in both cases. 

The basic sum-over states formula, Eq. (45), applies 
therefore to both OBCs and PBCs. In general, it is not 
very useful on practical grounds, since it would require 
the evaluation of slowly convergent sums. Nonetheless, 
Eq. (45) is instead essential to gather understanding into 
the physical meaning of rja^ , as will be shown below. 



CONDUCTIVITY 

Linear response 

So far, we have only discussed ground-state properties 
of our A''-electron system. Suppose now that it is subject 
to a small time-dependent perturbation contributing to 
the Hamiltonian the term: 



^ 1 - 

da;/(u;)-(Ae— * 



i^e'"*), (46) 



where A determines the "shape" of the perturbation and 
/ its amplitude. In order to get an Hcrmitian SH, we 
assume /(w) = f{—uj). We wish to measure the response 
to such perturbation by means of the expectation value 
of some observable B, i.e.: 



5{B) = - (*o|S|*o), 



(47) 



where ^ = + 5'^{t) is the perturbed time-evolved 
ground state. If we limit oiirselves to study terms which 
are linear in the response, it is enough to consider the 
single oscillatory perturbation: 



-iujt 



(48) 



whose response can be written, using the compact 
notations due to Zubarev [38-40], as: 



6{B) = ^{{{B\A))^ 



+ ((5|A))_,e^-*). (49) 



The quantity {{B\A))^ is by definition the linear response 
induced by the perturbation A at frequency w on 
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the expectation value (B). Straightforward first-order 
perturbation theory provides its exphcit expression as: 



(1'o|S|*„)(*«|i|'i'o) 



n#0 \ 

(^o|i|^„}(^„|^|^o} 
to + LOon + iri 



^ - Won + iri 



(50) 



where ujqi„ are the excitation frequencies of the 
unperturbed system, and the positive infinitesimal ry 
ensures causality. Expressions of the kind of Eq. (50) 
go under the name of Kubo formulae. 



Kubo formula for conductivity 

The conductivity tensor (7q,^(w) measures the current 
linearly induced by an electric field: — UafiSis- We 
therefore identify A with the potential of an electric field 
along (3, i.e. A = eSr^, and B with the current operator 
—eVa/L^- An important detail must be stressed at this 
point. The macroscopic field inside the; sample includes 
by definition screening effects due to the electronic 
system, while the perturbation 5H entering Eq. (29) — 
via the A operator — is the "bare", or unscreened one. 
This point will be discussed below (Sect. ); for the time 
being we simply identify screened and unscreened fields. 

The Kubo formula for conductivity is therefore 



(Ta/3(w) 



(51) 



this is correct within OBCs. but meaningless within 
PBCs, owing to the explicit presence of the position 
operator. This, however, makes no harm, since only its 
off-diagonal matrix elements are required: see Eq. (50). 
As usual, we may exploit the identity (\l/o|r|^'n) = 
i{^o\v\'^n)/^On- The Kubo formula becomes then 



o'afl(w) = lim —— 



(*o|^a|*n)(*n|^/3|*o) 



+ 



UJ + CJOn + if) 



W - iOOn + if] 



(52) 



We introduce a compact notation for the real and 
imaginary parts of the numerators in Eq. (52), i.e. 

^n,a/3 = Re (*o|Wa|«'™)('i'n|«/3|«'o), (53) 
In,a0 = lm{^o\Va\^n){^n\M^o), (54) 

which are symmetric and antisymmetric, respectively. 
Using then 



lim — - — = V— — inStx) 

r;->0+ X + iri X 



(55) 



and omitting the principal part, we separate for w > the 
symmetric and antisymmetric parts in the conductivity 
tensor as 



Reai+)(u;) = 



Re ay (a;) 



— d(w - Won) 



ire 

- ^ V' — 



On 



■ U) 



2 ■ 



(56) 



Sum rules 

At this point, we are ready to compare with the sum- 
over-states formulae for the tensor {rar^)c- In the present 
notations, we rewrite Eq. (45) as 



Re {rocr0)c 



-E 



2 

On 



Im {rar0)c = '^"'"^ 



(57) 



On 



A glance at Eq. (56) shows that 



Re {rar.h = ^ ^ ^ (58) 



Im {rccTf, 



ire^N 
2e^N 



Re 



,(-) 



(0). 



(59) 



Eq. (58) has been arrived at by SWM in 2000 [14], and 
Eq. (59) by Resta in 2005 [37]. First of aU, these identities 
show that {rar/s)c, defined here as a basic geometric 
feature, is indeed a measurable quantity (whenever it 
does not diverge). 

As emphasized throughout this work (rar^)c is a 
ground-state property, while the rhs of Eqs. (58) and 
(59) are properties of the system excitations, owing to 
the Kubo formula. Indeed, both Eqs. (58) and (59) 
look like the zero-temperature limit of a fluctuation- 
dissipation theorem, several forms of which are known in 
statistical physics [41, 42]: in the Ihs we have a ground- 
state fluctuation — see in particular Eq. (33) — while the 
ingredient of the rhs is conductivity (dissipation). 



Screened vs. unscreened field 

The Kubo formula for conductivity has been obtained 

identifying the A operator with Srjj, where £ is the 
macroscopic field inside the sample. In general this is 
not quite correct, since instead A = eSofp, where £q 
is the "bare" field, i.e. the field that would be present 
inside the sample in absence of screening. The latter 
originates from the two-body (electron-electron) terms 
in the potential V entering Schrodinger equation. The 
relationship between £ and Sq is not a bulk property. 
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and depends on the shape of the sample. Alternatively, 
it depends on the boundary conditions assumed for 
integrating Poisson equation: we refer to Ref. [12] for 
a thorough discussion. Whenever £ ^ £q, the sum rule 
in Eq. (58) must be modified. 

The ground-state fluctuations, as e.g. Eq. (33), 
are fluctuations of the macroscopic polarization, which 
in a finite sample induce a surface charge at the 
boundary. This in turn generates a depolarizing 
field, which counteracts polarization. Therefore the 
localization tensor (rar^)c depends on the sample shape, 
or equivalently on the boundary conditions assumed 
when taking the thermodynamic limit. The choice of 
PBCs in Eq. (29), however, implies £ = £o [43]; hence 
Eqs. (58) and (59) are correct as they stand within PBCs. 
This no longer holds within OBCs: in this case Eq. (58) 
needs to be modified, while Im (rar^)c = 0. 

Ideally the equality £ = £o corresponds to choosing 
a sample in the form of a slab, and to addressing the 
component of the fiuctuation tensor (rar/3)c parallel to 
the slab [12]; the thermodynamic limit amounts then to 
the infinite slab thickness. Owing to the long range of 
Coulomb interaction the order of the limits (first a slab, 
then its infinite thickness) is crucial. For instance, if the 
limit is taken instead by considering spherical clusters of 
increasing radius, the SWM fluctuation-dissipation sum 
rule, Eq. (58), assumes a different form: this is discussed 
in Refs. [36, 43]. The explicit form of the generalized 
sum rule is given therein. 

Last but not least, the effect leading to £ ^ £q 
within OBCs is a pure correlation effect. It originates 
from explicitly correlated wavefunctions, and does not 
occur within mean-field theories (Hartree-Fock and 
Kohn-Sham) [36, 43]. Within such theories, therefore, 
the sum rules hold in the simple form of Eqs. (58) 
and (59); the conductivity therein is the independent- 
particle conductivity ( "uncoupled" response in quantum- 
chemistry jargon). 



LOCALIZATION IN THE INSULATING STATE 

The basic tenet of the modern theory of the insulating 
state is that the localization tensor (rar^)c is the ground- 
state property which sharply discriminates — in the spirit 
of Kohn's seminal work [7, 8] — between insulators and 
metals. The real part of (rcir^)c remains finite in the 
thermodynamic limit in any insulator, while it diverges 
in any metal. 

The theory is very general, and has found applications 
to various different kinds of insulators: band insulators 
[44 47]; correlated (i.e. Mott) insulators, either by 
means of Hubbard-like model Hamiltonians [13, 48] or 
realistic ones [49]; "quantum Hall insulators" [37], and 
Anderson insulators [50]. As for Chcrn insulators and 
topological insulators, no explicit application of the 



present theory exists; nonetheless the work of Refs. 
[51, 52] implicitly shows that even in these cases the 
ground-state wavefunction is indeed localized in the sense 
of the present review. Most of these applications are 
reviewed in Sect. . 

The ultimate proof of the key property of Re {rarp)c is 
based on the SWM sum rule, Eq. (58). Since the tensor 
is real symmetric, it is enough to consider the diagonal 
elements (over its principal axes) 



Re {rara)c = 



hL 



3 /-cxj 



/ 

Jo 



dtu 



The /-sum rule yields 



Re craa(w). (60) 



f 

Jo 



doj Re cr„a(a;) = — 



8 2m^X3 ' 



(61) 



where cjp is the plasma frequency. Therefore the 
integral in Eq. (60) always converges at oo; its 
convergence/divergence is dominated by the small-w 
behavior of Re aaa (t^) • 

Suppose first that the spectrum is gapped, i.e. the 
spacing between the ground state and the first excited 
state stays finite in the thermodynamic limit. If the gap is 
Eg the conductivity vanishes for ui < Eg/H, and Eq. (61) 
yields 



duj 



< 



Re 



(T, 



Tre'^NE^ 
2meEg ' 



/•oo 

- / dujRe 
: Jo 



aa(w) 



(62) 



This inequality is due to SWM and clearly proves that 
Re {rarp)c is finite in any gapped insulator, as e.g. band 
insulators (considered in more detail in Sect. ). 

The main message of Kohn's 1964 paper, however, is 
that "insulating characteristics are a strict consequence 
of electronic localization (in an appropriate sense) and 
do not require an energy gap" . For any gaplcss material, 
the small-w behavior of Re <Jaa{'^) is the result of a 
competition between numerators and denominators in 
the Kubo formula, Eq. (52). Since we aim at a continuous 
function of w, the singularities in Eq. (56) must be 
smoothed: this can be done by keeping the "dissipation" 
rj finite while performing the thermodynamic limit 
first [53]. For a band metal the localization tensor 
diverges (see below). According to SWM, a gapless 
material is insulating whenever Re aaa{'-^) like a 
positive power of w, and metallic otherwise. The only 
example of gapless insulator considered so far is a model 
Anderson insulator in Id [50]. Simulations prove indeed 
that (a;^)c is finite therein (Sect. ). 
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Independent electrons 

For noninteracting electrons the potential V in 
Eq. (29) is the sum of identical one-body terms: 
V = J2iLi^{^i)- The many-electron Hamiltonian is 
separable and the exact ground state |^'o(*«)) is a Slater 
determinant of onc-particle orbitals (doubly occupied in 
the singlet case). At a mean- field level, the one-body 
potential V(r) includes electron-electron interaction in a 
selfconsistent way. In the Hartrcc-Fock (HF) framework 
the Slater determinant is regarded as an approximate 
many-electron wavefunction. Instead, in the density- 
functional framework the orbitals — called Kohn-Sham 
(KS) orbitals — are auxiliary quantities, individually 
devoid of physical meaning. In particular, their Slater 
determinant does not coincide with the many-electron 
wavefunction |^'o(k)), as a matter of principle. Therefore 
the exact localization tensor does not coincide, at least 
in principle, with the one obtained from the Slater 
determinant of KS orbitals. 

So much for the matters of principle. On practical 
grounds such difference is routinely disregarded (e.g. 
when dealing with polarization, magnetization [10-12], 
and more), given that it is not at all clear what is the 
relative importance of this "intrinsic" error, compared 
with the errors due to the choice of the functional 
itself. We therefore address here either the HF or the 
KS wavefunction |^'o(k)), having the form of a Slater 
determinant. 

Whenever the wavefunction is a Slater determinant, all 
ground-state properties can be explicitly cast in terms of 
the one-body density matrix 

N/2 

p(r,r') = 2P(r,r') = 2 ^ ^,(r)^*(r'), (63) 

where a singlet ground state is assumed, and ipj{r) are 
the occupied one-particle orbitals (either HF or KS); 
P(r, r') is the projector over the occupied manifold. 

The expression for the localization tensor is easily 
found within OBCs starting from Eq. (33) [36]: 

(r„r^)c = ^ Jdrdr' (r - r')a(r - r')^ |P(r, r')p. (64) 
If we define the complementary projector 

g(r,r') = (5(r-r')-P(r,r'), (65) 
an equivalent expression is [44] 

(rar/3)c = |:Tr {r^PrpQ}, (66) 

where "Tr" is the trace over the single-particle Hilbert 
space (not on Cartesian indices). 



If we consider a cluster, cut out of a crystalline solid, 
Eq. (64) becomes in the large- limit 

{vaV^U = ^ [ dr [ dr' (r - r'Ur - r')^ |P(r,r')|^ 

c J cell J all space 

(67) 

where is the number of electrons per crystal cell. 
According to the discussion in Sec. , we need not to 
worry about shape issues in taking the limit; we also 
notice that the density matrix, Eq. (63), is independent 
of the boundary conditions (either OBCs or PRCs) in the 
large- A^ limit. 

Band insulators and band metals 

As observed, Eq. (67) holds for a crystalline solids. 

Therefore the inner integral on the rhs must converge 
for a band insulator, and must diverge for a band 
metal. This is confirmed by the well known fact 
that the asymptotic behavior of P is qualitatively 
different in insulators and in metals. In the former 
materials, in fact, P(r, r') decays exponentially [54-57] 
for large values of r — r': therefore the integral converges 
and the localization tensor is finite. In conducting 
materials, instead, P(r, r') decays only polynominially, 
and the inner integral diverges. This divergence can 
be explicitly verified for the simplest conductor of all, 
namely, the noninteracting electron gas, whose density 
matrix is exactly known in analytic form [44, 58]. 
Therefore the localization tensor, when expressed in the 
form of Eq. (67), measures in a perspicuous way the 
"nearsightedness" [59] of the electron distribution. Such 
measure is qualitatively different in insulators and in 
metals. 

The onc-particle orbitals (cither HF or KS) in a 
crystalline solid have the Bloch form. We therefore may 
wish to replace the orbitals in the expression for P(r, r'), 
Eq. (63), 

<^i(r) ^ Vnq(r) = e^'»--u„q(r), (68) 

where n is the band index and q is the Bloch vector. We 

stress that PBCs are at the very root of Bloch theorem. 
If the orbitals are normalized to one over the crystal cell 
of volume V^, the ground-state projector in insulating 
crystals is 

V ""^'^ r 

P(r,r') = / dqV„q(r)Cq(r') (69) 

i^^r ^1 Jbz 

V ""^'^ r 

= TT^E / c^qe^'>("^-"-)u„q(r)<<,(r'), 

where Nc/2 is the number of occupied bands, and the 
integral is taken over the Brillouin zone (equivalently, 
over the reciprocal cell). 
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Using Eq. (69), the localization tensor in a band 
insulator becomes (for double occupancy) 



dimensionless and simply related to the Chern number, 
defined as: 



(r r.) - / da (s^l^u \—u ) C, = -- f dq J] Im ( A^„q| (73) 

\ Tl 



-^{Un^\-^Un'ci){-^Un'ci\Un<i) \ -(70) 
n,n' oi P I 

The proof is given in Refs. [15, 44], and will not be 
repeated here. 



Wannier functions 

Expressions similar to Eq. (70) also enter the Marzari- 
Vanderbilt theory of maximally localized Wannier 
functions [60]; the relationship is 



(the formula here is given for single band occupancy). 

Band insulators where the Chern invariant is nonzero 
are called "Chern insulators" (normal insulators other- 
wise). The possiible existence of Chern insulators has 
been pointed out by Haldane in 1988 [61], by means of 
a remarkable model Hamiltonian in 2d, much studied af- 
terwards [51, 62, 63], and illustrated below (Sect. ). We 
emphasize that the nonvanishing of the Chern number 
prevents the existence of exponentially localized Wan- 
nier functions [51, 64], at variance with normal insulators 
where the "maximally localized" Wannier functions [60] 
are exponentially localized [65]. 



. ^ 2 



(71) 



Here Vti indicates the gauge-invariant part of the 
quadratic spread of the Wannier functions, as in Ref. 
[60] and in the subsequent literature. 

The trace in Eq. (71) is a lower bound (and not a 
minimum in dimension d > 1) for the spherical second 
(cumulant) moment — a.k.a. quadratic spread — of the 
Wannier functions, averaged over the sample. 

We stress that Eqs. (69) and (70) make sense only 
insofar the Fermi level falls in a gap, in which case (rQr^)c 
is always finite. If we vary the Hamiltonian continuously, 
allowing the gap to close, then {rarp)c diverges; the 
quadratic spread of the Wannier functions diverges as 
well [36]. 



Chern insulators 

Within OBCs the localization tensor is always real: 
this is perspicuous in Eqs. (64) and (67). Instead the 
PBC expression of Eq. (70) is naturally endowed with an 
imaginary part: 



Im {rarfi)c 



2K 



(27r)37Vc 



d 
dqa 



dqi3 



Wnq).(72) 



A necessary condition for this to be nonzero is the 
absence of time-reversal symmetry. 

The second line of Eq. (72) has (in 3d) the dimensions 
of an inverse length, and is quantized in units of tt times a 
reciprocal vector. We call this reciprocal vector "Chern 
invariant"; in fact it is a Chern number (Sect. ) in q- 
space and in appropriate units. In 2d the BZ integral is 



LOCALIZATION IN DIFFERENT KINDS OF 
INSULATORS 

Small molecules 

The modern theory of the insulating state clearly 
addresses extended systems, i.e. the TV — >• oo limit; 
indeed it makes little sense to ask whether a small 
molecule is insulating or conducting. Nonetheless the 
concepts of localized/delocalized electronic states is of 
the utmost importance in quantum chemistry as well, 
notably in relationship to aromaticity. 

The tensor {rarp)c within OBCs is always real 
symmetric. If the ground-state wavefunction is a Slater 
determinant, then the trace of the tensor at finite N has 
the meaning of a lower bound for the quadratic spread of 
the Boys localized orbitals, averaged over all the occupied 
orbitals [36]. 

The small-A/' version of the main concepts of the 
present review (in their OBCs fiavour [36]) has been 
adopted in quantum chemistry by Angyan [66, 67]. 
Besides providing HP calculations of {roir0)c for a sample 
of small molecules, Angyan even provides experimental 
values drawn from compilations of the dipole oscillation- 
strength distributions: basically, from Eq. (58). 



Band insulators 

The theory warrants that the localization tensor 
is finite in any insulator. However, quantitative 
calculations for both model tight-binding Hamiltonians 
and realistic solids within density-functional theory have 
been used to illustrate the theory and to identify trends. 
For instance one expects much smaller diagonal elements 
{'''af'a)c for strong (i.e. large-gap) insulators than for 
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weak (small-gap) ones. This is also suggested by the 
SWM inequality, Eq. (62). 

Let us start with a simple tight-binding (a.k.a. Hiickel) 
Hamiltonian in Id: 

^ = El (-1)'^ 4^^- - + H-c.) ] (74) 

jo- 

where t > is the first neighbor hopping {(3 = —t in 
most chemistry literature) and H.c. stays for Hermitian 
conjugate. This toy model schematizes a binary ionic 
crystal; the band structure is 



e{q) = ±VA2 +4*2 cos2 qa/2, 



(75) 



where a is the lattice constant and q is the Bloch vector. 
The gap is equal to 2A; at half filling the system is always 
insulating except for A = 0. The squared localization 
length (within OBCs) is the tight-binding version of 
Eq. (64), i.e. 



a 

AN 



N 



(76) 



This is a monothonical function of A; it is easily verified 
that it vanishes in the extreme ionic case {t = 0). In the 
metallic case (A = 0) the ground-state projector has a 
simple analytical form: 



Pjj = 2' ^ ^ ^^^^ 1-^' ~ ■^l = 



P / = 



(-ir 



~ n{2s + l) l/-j| = 25 + l, (77) 

which clearly implies divergence of Eq. (76). At any finite 
within OBCs Eq. (76) leads to a finite {x'^)c value; 
however Eq. (77) suggests that in the metallic case {x'^)c 
diverges linearly with N. This has been verified by actiial 
simulations, even when A ^ but the Fermi level is not 
in the gap [50]. 

Other simulations [47] have addressed dimerized 
chains, i.e. A = but alternant hoppings in Eq. (74). 
While nothing relevant occurs within PBCs, partly filled 
end states within OBCs at some fillings are at the root 
of some noticeable features. 

The first ab-initio study (in 2001) addressed several 
elemental and binary cubic semiconductors at the KS 
level [44] . The tensor is real and isotropic. The computed 
(a;2)c (Fig. 5) is smaller than 3 bohr^ in all the materials 
studied: the ground many-body wavcfunction is therefore 
very localized in this class of materials. The SWM 
inequality was also checked, and found to be well verified 
using both the theoretical KS gap and the experimental 
one (the latter is typically larger). 

Other studies have addressed the ferroelectric per- 
ovskites in their different (cubic and noncubic) structures 
[45], and some model Hamiltonians in Id and 2d [46]. 



o 



• egcpeiinn^gap 



as* .♦^a» 



ft 1 2 S 4 S 6 7 



9 10 11 12 W 1S 1£ 



FIG. 5: Diagonal element of the KS localization tensor vs. the 
inverse direct gap (theoretical and experimental), for several 
elemental and binary semiconductors (from Ref. [44]) The 
points corresponding to Si and Ge with the theoretical gaps 
are out of scale. Ftom Ref. [44]. 



Correlated (Mott) insulators 

Starting from the noninteracting Hamiltonian of 
Eq. (74) and augmenting it with an on-site repulsive term 
we get the two-band Hubbard model 



=E 

jo- 



[(-l)^Act c,.-t(ct c,+i.+H.c.)]- 



j 



(78) 

The explicitly correlated ground-state wavcfunction has 

been found by exact diagonalization [13], and the 
corresponding {x'^)c has been computed as a function 
of U for fixed t/A = 1.75. The results are shown 
in Fig. 6 in dimcnsionlcss units; it turns out that 
there is only one singular point U = 2.27t, where {x^)c 
diverges. Indeed, it has been verified that at such 
value the ground-state becomes degenerate with the first 
excited singlet state, i.e. the system is metallic. The 
singular point is the fingerprint of a quantum phase 
transition: on the left we have a band-like insulator, and 
on the right a Mott-like insulator. The two insulating 
states are qualitatively different; by adopting the modern 
jargon, nowadays we could say that they are topologically 
distinct. The static ionic charges (on anion and cation) 
are continuous across the transition, while the dynamical 
(Born) effective charge on a given site changes sign [68]. 
Other studies of the localization tensor within the same 
Hubbard model can be found in Ref. [69] . 

The transition from a band metal to a Mott insulator 
has been studied in a model linear chain of Li atoms 
by Vetere et al. [49]. At a mean- field level the infinite 
chain is obviously metallic at any lattice constant a, since 
there is one valence electron per cell. However the mean- 
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FIG. 6: Squared localization length for the Hamiltonian in 
Eq. (78) at half filling for t/A= 1.75. The system undergoes 
a quantum phase transition from band-like insulator to Mott- 
like insulator at U/t = 2.27. Prom Ref. [13] 



field description becomes inadequate at large a, where 

the electrons localize and the system becomes a Mott 
insulator. If electron correlation is properly accounted for 
at any a, the system undergoes a sharp metal-insulator 
transition at a critical a. 

The calculations addressed linear Lijv systems {N up 
to 8), where the finite size prevents a sharp transition; 
the tradeoff is that full configuration interaction was 
affordable with 6 atomic orbitals per site (yielding 
more than 10^ symmetry-adapted Slater determinants). 
The wavefunction of Vetere et al. is therefore exempt 
from any bias insofar as the treatment of correlation is 
concerned, although its quality is determined by the basis 
set. A study of the longitudinal component (a;^)c of the 
localization tensor indicates rather clearly the occurrence 
of the metal-insulator transition at a ~ 7 bohr; other 
indicators give concordant results [49]. For comparison, 
the nearest-neighbour distance in 3d metallic lithium is 
5.73 bohr. 



Disordered (Anderson) insulators 

We start from the same Hamiltonian as in Eq. (74), and 
we replace the ordered string (—I)-' by a random string 
of ±1, chosen with equal (and uncorrclated) probability. 
This system models a random binary alloy at 50% 
concentration. It is well known both from analytical 
arguments and actual simulations that its spectrum is 
gapless [6, 70]. The density of states for both the ordered 
and disordered systems are shown in Fig. 7, and confirm 
the expected features. The band structure of Eq. (75) 
yields obviously a gapped density of states; at the band 
edges it shows van Hove singularities, which in Id have 
the character of l/\/e divergences. As discussed above, 
the system is insulating at half filling and conducting 
otherwise. The disordered system, instead, is gapless 
and nonetheless insulating at any filling. In fact, this 
model Hamiltonian describes a paradigmatic Anderson 
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FIG. 7: Density of states (arbitrary units) for a model binary 
alloy in Id. The crystalline (band) case corresponds to the 
Hamiltonian of Eq. (74) with A = 0.25 and t = 1. The 
disordered (Anderson) case corresponds to a random choice 
of the anion/cation distribution. 



insulator in Id. 

The conventional theory of transport focusses on the 
nature of the one-particle orbitals at the Fermi level; in 
Anderson insulators these are localized, thus forbidding 
steady state currents [4]. More than fifty years of 
literature have been devoted to investigate Anderson 
insulators under the most diverse aspects. [6, 70-72]. 

At variance with such wisdom, a recent work has 
addressed this paradigmatic Anderson insulator from the 
nonconventional viewpoint of the modern theory of the 
insulating state [50]. In the spirit of Kohn's theory 
the invidual Hamiltonian eigenstates become apparently 
irrelevant, while the focus is on the many-electron ground 
state as a whole. The squared localization length (a;^)c 
has been computed within OBCs from Eq. (76), and 
found to be finite, as expected. Nonetheless its value is 
about 20 times larger than the one for the band insulator, 
at the same value of the parameters (i.e. A = 0.25, t = 
1). This reflects the fact that the scattering mechanisms 
are profoundly different: incoherent (Anderson) versus 
coherent (band). In the latter case, the Hamiltonian 
eigenstates are individually conducting but "locked" by 
the Pauli principle if the Fermi level lies in the gap. 



Quantum Hall insulators 

The results of a typical quantum Hall (QH) experiment 
are shown in Fig. 8. The most perspicuous feature is 
that whenever the 2d electron fluid is in the QH regime: 
(i) the transverse resistivity is quantized; and (ii) the 
longitudinal resistivity vanishes. While the plots are 
about the integer QH effect, the same two features are 
common to the fractional QH effect as well. 
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FIG. 8: The results of a typical quantum Hall experiment. 



We switch from resistivity to conductivity; in isotropic 
2d-systems the relationship is 



Pii 



0-11 



Pl2 



0'12 



(79) 



hence in he QH regime an — and <7i2 = — l//3i2 
is quantized. Since the static longitudinal conductivity 
vanishes, the system is by definition insulating: this is 
why we speak of "QH insulators" . 

We are going to show that electron localization — as 
defined in the modern theory of the insulating state — 
is the common cause for both the vanishing of the dc 
longitudinal conductivity and the quantization of the 
transverse one. Therefore to predict whether the dc 
transverse conductivity is quantized, it is enough to 
inspect electron localization in the ground state. This 
outstanding result was proved in 2005 [37]. 

Our system is a square of size LxL, where the presence 
of a macroscopic B field forbids simple PBCs; this, 
however, is no serious problem. We choose the Landau 
gauge for the magnetic vector potential in Eq. (29), 
and we impose the magnetic boundary conditions which 
are customary in the QH literature [73]; this requires 
the total flux BL^ to be an integer times $0) where 
$0 = hc/e is the flux quantum. 

When switching from 3d to 2d, many of the above 
formulae require obvious modifications. For an isotropic 
system, the 2d analogue of Eqs. (58) and (59) can be 
written as 



(80) 



Re 0-12 (0) = 



2e^N 



Im {xy),. (81) 



A glance at the sum-over-states expression for {rarp)c, 
Eq. (45), shows that whenever Re {rarj3)c converges. 



then Im {rarfj)^ converges as well. However, we are going 
to show that the latter may only converge to quantized 
values in the present 2d case. 
By definition we have 



Im {xy)c = ^Im (9i*o(0)|%*o(0)). 

In the large- L limit, wc may write 



(82) 



Im (9i«'o(0)|a2*o(0)) 

2 

Im / dKi 



= (- 

\2tt 







(83) 

^dK2 (ai*o(«)|92*o(K)). 



The dimensionless integral is over a closed surface (a 
torus), given the magnetic boundary conditions. It is 
therefore equal to — ttCi, where Ci is the Chern number 
of the first class (Sect. ); hence 

and the transverse conductivity, Eq. (81) becomes 

We have thus arrived at the famous Niu-Thouless- 

Wu formula [73] which holds for both the integer 
and fractional QH effect. The original derivation was 
based on an analysis of the Green function, under the 
hypothesis that the system has a Fermi gap; in the 
present derivation the presence of a Fermi gap becomes 
apparently irrelevant, since the localization tensor is 
a pure ground-state property. Underlying the modern 
theory of the insulating state, however, is a fluctuation- 
dissipation theorem, relating the ground state to the 
excitations of the system. 

The main message of Ref. [37] (and of the 
present Section) can be summarized by saying that 
any 2d insulator displays a quantized transverse 
conductance (nonvanishing only in absence of time- 
reversal symmetry). 

An electron fluid is kept in the QH regime by disorder, 
and analytical implementations of the present formulae 
are obviously not possible. Nonetheless, in order to 
illustrate how the theory works, it is expedient to 
consider the academic case of noninteracting electrons 
in a flat substrate potential. If the first Landau level 
is fully occupied (i.e. at filling 1^=1) the electron 
density is uniform and equal to no = l/(27ri!^), where 
i = {hc/eBy/^ is the magnetic length; the modulus of 
the density matrix has the simple expression (for single 
occupancy) 



|p(r,r')| = |P(r,r')|=noe 



-|r-r'|V(«=^) 



(86) 



and the cumulant second moment is clearly finite. The 
trace {x^)c + {y'^)c of the locahzation tensor is in fact 
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equal to , the squared magnetic length. If the 
filling is fractional the density cannot be uniform (for 
noninteracting electrons), and, more important, the 
cumulant second moment diverges (see Ref. [37] for 
details): the system is not insulating, and its transverse 
conductivity is not quantized. 

Chern insulators 

The QH state of matter, discovered in 1980 and 
addressed in the previous Section, provides the first 
example of a quantum state which is topologically 
distinct from all previously known states of matter. 
This is the reason why its macroscopic quantization 
properties are very robust ("topologically protected"), 
and insensitive to small changes in material parameters. 
The debut of geometrical concepts in electronic structure 
theory dates from two milestone papers: TKNN in 
1982 [24] and Niu-Thouless-Wu in 1985 [73]. These 
papers addressed a 2d electron fluid (noninteracting and 
interacting, respectively) in the presence of a macroscopic 
magnetic field: the Hamiltonian, therefore, cannot be 
lattice periodical. 

A subsequent breakthrough on the theory side is the 
Haldane model Hamiltonian [61]: this can be considered 
as the precursor of modern topological insulators [31]. Its 
trademark is quantized Hall conductance in absence of a 
macroscopic magnetic field. 

The model is comprised of a 2d honeycomb lattice 
with two tight-binding sites per primitive cell with site 
energies ±A, real first-neighbor hopping ti, and complex 
second- neighbor hopping t^e^^'^ ^ as shown in Fig. 9. 
Within this two-band model, one deals with insulators 
by taking the lowest band as occupied. The appeal of the 
model is that the vector potential and the Hamiltonian 
are lattice periodical and the single-particle orbitals 
always have the usual Bloch form. Essentially, the 
microscopic magnetic field can be thought as staggered 
(i.e. up and down in different regions of the cell), but its 
cell average vanishes. 

For a 2d lattice-periodical Hamiltonian the Chern 
number C\ has been defined before, Eq. (73). As a 
function of the flux parameter this system undergoes 
a transition from a normal insulator {C\ = 0) to Chern 
insulator (|Ci| = 1). 

Like in the QH case discussed in Sect. : (i) the 
imaginary part of the localization tensor is related to 
the Chern number by Eq. (84); and (ii) the trace of the 
localization tensor is finite. This is confirmed by the 
simulations of Ref. [51], where the actual value of fli, 
Eq. (71), is computed. We remind that, despite fli being 
finite, localized Wannier functions (with finite quadratic 
spread) do not exist in Chern insulators. 

The Haldane model insulator and the QH insulator 
both display quantized transverse conductivity, and both 
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FIG. 9: Four unit cells of the Haldane model [61]. Filled 
(open) circles denote sites with iJo = —A (-I-A). Solid 

lines connecting nearest neighbors indicate a real hopping 
amplitude t\\ dashed arrows pointing to a second- neighbor 
site indicates a complex hopping amplitude tie^'^ . Arrows 
indicate sign of the phase 4> for second-neighbor hopping. 

are localized in the sense of the modern theory of the 
insulating state. It is worth pointing out, though, that 
the decay of the density matrix is qualitatively different: 
exponential in the Haldane case [51] vs. Gaussian in the 
(noninteracting) QH case, Eq. (86). 

No microscopic realization of a Chern insulator (in 
absence of a macroscopic magnetic field) is known to 
date. A mesoscopic Chern insulator, in the same spirit as 
the Haldane model, was synthesized in 2008 [74]. Chern 
insulators remained a curiosity of academic interest only 
for many years. In 2005 it was realized that the Haldane 
model is the forebear of a completely new class of 
insulators, called "topological insulators" [25, 31]. 

Topological insulators 

Although Chern insulators do not exist in nature, 
a different kind of topological insiilators does exist in 
nature in both 2d and 3d [26 31]; their macroscopic 
quantization properties are topologically protected, and 
insensitive to small changes in material parameters, 
exactly like in the QH case. 

Modern topological insulators are characterized by 
novel invariants called Z2 (close relatives of the Chern 
number in 2d and of the Chern invariant in 3d). The 
paradigmatic Z2 topological insulator is the Kane-Mele 
model Hamiltonian in 2d [25]. Most important, at 
variance with Chern insulators, the modern topological 
insulators are time-reversal symmetric. 

Since spin-orbit interaction plays a major role, the 
definition of localization tensor needs to be augmented 
to explicitly include the spin coordinates; this seems 
to be rather straightforward, and is deferred to further 
publications. After this is done, the scope of the modern 
theory of the insulating state would include topological 
insulators. 

Obviously, no investigation of this kind exists yet. 
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Nonetheless, a closely related investigation concerning 
the Wannier functions in a Z2 insulator is appearing these 
days [52] . The results confirm that the ground state of 
a topological insulator is localized in the sense of the 
modern theory of the insulating state: ^li, Eq. (71), is 
in fact finite. Furthermore — at variance with the Chern- 
insulator case — localized Wannier functions (with finite 
quadratic spread) do exist, although a generalization of 
the Marzari-Vanderbilt theory [52, 60] is needed. 



CONCLUSIONS 

According to the milestone 1964 Kohn's paper [7], 
the organization of the electrons in their ground state 
is qualitatively different in insulators and metals. This 
paper provided the foundations of the "theory of the 
insulating state" . The modern developments of the 
theory started many years later (1999) with the work 
of Resta and Sorella [13], and continue to these days. 

Here we provide a comprehensive presentation of 
the modern theory. In so doing, we emphasize the 
geometrical properties of the many-body wavefunction, 
which were somewhat hidden in most of the original 
liter atiirc. 

We focus on the linear response of the many-body 
wavefunction to an infinitesimal "twist" (or "flux") in 
the electronic Hamiltonian. This is the ground-state 
property which discriminates insulators from metals. 
The distance between the (inflnitesimally) twisted and 
untwisted ground states, measured by the quantum 
metric tensor in the thermodynamic limit, is finite in 
insulators and divergent in metals. 

Technically, in our formulation the twist has the 
dimensions of an inverse length; the quantum metric 
tensor has therefore the dimensions of a squared length 
and is extensive. The key intensive quantity is the 
metric tensor per electron; this is generally called 
"localization tensor" or "second cumulant moment" of 
the electron distribution. It is a global property of 
the electronic ground state, apparently unrelated to the 
system excitations. In some cases (discussed above) 
the localization tensor is naturally endowed with an 
imaginary antisymmetric part, whose physical meaning 
is also outstanding. 

The modern theory of the insulating state encompasses 
all kinds of known insulators. Here we discuss, using 
uniform concepts and notations, the most diverse cases: 
open and periodic boundary conditions; intciracting 
and noninteracting electrons, crystalline and disordered 
systems; time-reversal symmetry present and absent. We 
even briefly comment on topological insulators. 

Work partly supported by the ONR grant N00014-07- 
1-1095. 
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